OKHEP-99-05 

Solution of Schwinger-Dyson Equations for PT-Symmetric 

Quantum Field Theory 

Carl M. Bender/,* Kimball A. Milton,2t and Van M. Savage^^ 
^Department of Physics, Washington University, St. Louis, Missouri 63130, USA 
'^Department of Physics and Astronomy, University of Oklahoma, Norman, Oklahoma 73019, 



In recent papers it has been observed that non-Hermitian Hamiltonians, 
such as those describing igcp^ and —g4>'^ field theories, still possess real pos- 
itive spectra so long as the weaker condition of VT symmetry holds. This 
allows for the possibility of new kinds of quantum field theories that have 
strange and quite unexpected properties. In this paper a technique based 
on truncating the Schwinger-Dyson equations is presented for renormalizing 
and solving such field theories. Using this technique it is argued that a —gcp'^ 
scalar quantum field theory in four-dimensional space-time is renormalizable, 
is asymptotically free, has a nonzero value of (0|(/>|0), and has a positive defi- 
nite spectrum. Such a theory might be useful in describing the Higgs boson. 



USA 

(February 1, 2008) 



Abstract 



11.15.Pg, ll.30.Qc, 25.75.-q, 3.65.-W 



Typeset using REVT^ 



*electronic mail: 



cmb@howdy.wustl.edu 



^electronic mail: 



milton@mail.nhn.ou.edu 



■^electronic mail: 



vmsavage@hbar.wustl.edu 



1 



] 



I. INTRODUCTION 

It has recently been observed that quantum mechanical theories whose Hamiltonians 
are PT-symmetric have positive definite spectra even if the Hamiltonian is not Hermitian. 
A class of such theories that has been studied extensively is defined by the Hamiltonian 

H = p^-{ixf iN>2). (1.1) 

It is believed that the reality and positivity of the spectra are a direct consequence of VT 
symmetry. 

The positivity of the spectra for all is an extremely surprising result; it is not at 
all obvious, for example, that the Hamiltonian H = p"^ — corresponding to = 4 has 
a positive real spectrum. To understand this result it is necessary to define properly the 
boundary conditions in the corresponding Schrodinger equation. 

For the Hamiltonian in Eq. ( |1 . 1| ) the Schrodinger differential equation corresponding to 
the eigenvalue problem Hip = Eip is 

- i)"{x) - {ix)^ip{x) = EiP{x). (1.2) 

The boundary conditions for this equation are discussed in detail in Ref. [Q. There, it was 
shown how to continue analytically in the parameter A^ away from the harmonic oscillator 
value N = 2. This analytic continuation defines the boundary conditions in the complex- a; 
plane. The regions in the cut complex-x plane in which iplx) vanishes exponentially as |a;| — >■ 
oo are wedges. In Ref. [|l| the wedges for A^ > 2 were chosen to be analytic continuations of 
the wedges for the harmonic oscillator, which are centered about the negative and positive 
real axes and have angular opening |. For arbitrary N > 2 the anti-Stokes' lines at the 
centers of the left and right wedges lie below the real axis at the angles 
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/A^-2\ TT 

VA^ + 2/ 2' 

^nght - - [j;^) 2 ■ (1-3) 

The opening angle of these wedges is -^^^ In Ref. IH the time-independent Schrodinger 
equation was integrated numerically inside the wedges to determine the eigenvalues to high 
precision. 

The quantum mechanical Hamiltonian in Eq. ( [I.ID has additional remarkable properties. 
For example, for all A^ > 2 the expectation value (0|x|0) of the position operator x in the 
ground state is nonzero. This surprising result is true even when A^ = 4 |Q. 

These results for quantum mechanics raise some interesting questions regarding quantum 
field theory. In particular, does the self-interacting scalar quantum field theory defined by 
the Lagrangian 

£ = 1(90)2 + lmV-^W^ (1.4) 
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have a positive definite spectrum and a nonvanishing value of (O|0|O) for all N > 21 We 
believe that the answer to this question is yes. Because of these properties, we believe that 
when = 4 the resulting quantum field theory in four- dimensional space-time could serve 
as a good description of the Higgs particle. As we argue in this paper, the —gcj)'^ theory is 
particularly advantageous because, like the conventional gcf)'^ it has a dimensionless coupling 
constant, but unlike the conventional theory, it is asymptotically free and is thus not a trivial 
theory. 

The question of how to determine the properties of PT-symmetric non-Hermitian quan- 
tum field theories has already been examined. As we will see in this paper, conventional 
Feynman diagrammatic perturbation theory is not adequate for studying these theories. 
Thus, in previous work the perturbative approach that was used was to take N = 2 + 5, 
where 5 is treated as a small parameter. While some interesting results regarding parity 
violation and supersymmetry were obtained, unfortunately this perturbative scheme 
has a severe drawback: It is not known how to carry out the renormalization procedure 
required to understand higher-dimensional theories. 

Why is it that Feynman diagrams cannot be used to perform calculations in theories 
such as —g(t)^l As we have already stated, in this theory (O|0|O) ^ 0. There is no way to 
obtain this result using the standard Feynman rules; one cannot obtain a one-point Green's 
function using four-point vertex amplitudes. Indeed, as we will show in Sec. |T| in the context 
of zero-dimensional theories, the standard Feynman rules are incorrect for this theory. 

In this paper we perform a systematic truncation of the Schwinger-Dyson equations as 
a calculational procedure. This idea has already been applied in a simple context to obtain 
Green's functions and energy levels in conventional quantum mechanical problems . Trun- 
cating the Schwinger-Dyson equations is an inherently variational approach; including more 
and more of the higher Green's functions is equivalent to enlarging the space of variational 
parameters. In a recent study of the PT-symmetric, non-Hermitian quantum mechanical 
Hamiltonian in Eq. ( |1.1D variational methods were found to be extremely accurate [Q. 

This paper is organized as follows: In Sec. |I| we review the general approach used in 
this paper. We show how to derive the Schwinger-Dyson equations using simple functional 
methods and we explain our truncation procedure. In Sec. |T| we examine the numerical 
accuracy of our truncation method in the context of zero-dimensional field theory. We study 
the massless zero-dimensional version of Eq. ( |1.4| ) in very high order in the truncation process 
for the cases = 3 and = 4. In the case of a zero-dimensional massive theory we show 
that the Feynman rules are inapplicable. In Sec. |V| we examine the theory in Eq. ( |1.4D in 
one-dimensional space-time (quantum mechanics). Finally, in Sec. 0we apply our methods 
to field theory in arbitrary dimension. As an application of our procedure we calculate the 
Callan-Symanzik function (3{g) for a four-dimensional —gcj)'^ theory for small g. We show 
that to leading order (3{g) is negative and thus the theory is asymptotically free. This result 
is in distinct contrast with the result for a conventional gcj)'^ theory. 

II. ELEMENTARY DERIVATION OF THE SCHWINGER-DYSON EQUATIONS 

The Schwinger-Dyson equations are an infinite set of coupled equations relating the 
Green's functions of a quantum field theory. In this section we derive the Schwinger-Dyson 
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equations using elementary formal functional methods. 

We begin with Eq. (lA) and append a term that represents the coupling of the field 0(a:) 
to an external c-number source J{x): 

^ = l{d<Pr + IrnY - j^{^'^>r - (2.1) 

This Lagrangian represents a self-interacting scalar quantum field theory in D-dimensional 
Euclidean space-time. If we vary the action with respect to (\)[x) we obtain the field equation: 

- (9V(x) + m^(\){x) - ig{v\){x)f-^ = J{x). (2.2) 

Next, leaving the source turned on, we take the expectation value of the field equation 
in the vacuum state of the theory |0) and divide by the vacuum- vacuum functional 
Z[J] = (0|0): 

- a^G'/\.) + m-G'A.) - 9'" '°'^ (01(1""°' = J(^), (2.3) 
where G['^\x) is the one-point Green's function in the presence of the external source: 

e,.,,) . M. (2.4) 



Note that the function J{x) appears alone on the right side of Eq. ( P73| ) because it is a 
c-number and therefore can be factored out of matrix elements. 

The objective is now to use Eq. ( |2.3|) to calculate the Green's functions of the theory. 
The (connected) Green's functions in the presence of the source J are defined as functional 
derivatives of the logarithm of Z[J] with respect to the source J{x): 

Gi'\x„ X2, . . . , x„) = —— HZ[J]). (2.5) 

dJ[Xi)dJ{X2) ■ ■ ■Oj{Xn) 

To obtain the standard connected Green's functions of the theory (the connected part of the 
vacuum expectation value of the time-ordered product of the fields; that is, the sum of the 
connected n-point Feynman diagrams) we then turn off the source: 

Gn{Xi, X2,. .. , Xn) = G^^^Xi 

) r-. ■ (2-6) 

Note that turning off the source restores translation invariance. As a result, the one-point 
Green's function Gi is a constant independent of x, the two-point Green's function depends 
only on the difference x — y, G2{x, y) = G2{x — y), G^i^x, y, z) = G^^x — y,x — z) depends 
on two differences, and so on. 

Before one can proceed, one must express the third term in Eq. ( ^.3| ) in terms of the con- 
nected Green's functions of the theory. To do this we recall that functionally differentiating 
with respect to J{x) is equivalent to inserting 0(a;) in matrix elements [^]. 

Let us consider the simple case iV = 3. In this case it is necessary to calculate the 
quantity (O|0^(a;)|O). To do so we begin with Eq. ( p.4| ) multiplied by Z[J]: 
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G[^\x)Z[J] = (0|(/)(x)|0). (2.7) 

Taking the functional derivative of this equation with respect to J{x) gives 

[G[''\x)]^Z[J] + G^^\x,x)Z[J] = {0\(f)\x)\0). (2.8) 

Hence, we can ehminate (O|0^(a;)|O) from Eq. ( p.3| ) to obtain 

-d^G[''\x)+m'G[^\x)+gt{[G[''\x)]^ + Gi^\x,x)) = J{x). (2.9) 

We now obtain the first of the Schwinger-Dyson equations by setting J = (turning off the 
source) : 

m'^Gi + gi [gI + 6*2(0)] = 0. (2.10) 

Remember that by translation invariance Gi is a constant, so that its derivative vanishes 
and that 6*2(0) = G2{x — x) = G2{x,x). 

To obtain the second of the Schwinger-Dyson equations for = 3 we take a functional 
derivative of Eq. (p.9|) with respect to J{y), 

- d'Gi''\x,y)+m'Gi''\x,y)+gt [2G[''\x)Gi'\x,y) + Gi'\x, x,y)] = 5{x - y), (2.11) 

and then set J = in this equation: 

- 9^G2(x -y)+ n?G2{x - y) + gi [2GiG2ix - y) + 6*3(0, x - y)] = 5{x - y). (2.12) 

If we continue the process of functionally differentiating with respect to J and setting 
J = 0, we obtain the infinite tower of coupled differential equations known as the Schwinger- 
Dyson equations. For example, the third in the sequence is 

—d'^G-si^x — y,x — z) + m^G^{x — y,x — z) 

+gi[2GiG3{x -y,x-z) + 262(0; - ^)6'2(x - y) + Gi{0, x-y,x- z)]=0. (2.13) 

Note that these Schwinger-Dyson equations are incomplete in the sense that there are 
too many unknowns. The first equation contains Gi and G2, the second contains 6*1, 6*2, 
and 6*3, and so on. Thus, each new equation contains a new unknown Green's function 
and the system never closes. However, we can force the system to close by truncating the 
sequence of coupled equations and setting Gn+i = in the last equation. We will use this 
method throughout the remainder of the paper. 

As a second example we derive the first four Schwinger-Dyson Equations from Eq. ( [2.3[ ) 
for the case = 4. Using the same approach as we did for the case A^ = 3, we begin by 
reexpressing {0\(j)^{x)\0). We do this by taking the functional derivative of Eq. ( |2.8| ) with 
respect to J{x) to obtain 

[G[^\x)fZ[J]+3G[^\x)Gi^\x,x)Z[J]+G^3^\x,x,x)Z[J] = (O|0'^(x)|O). (2.14) 

Substituting this result into the A^ = 4 version of Eq. (|2.3|), we have 
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- d^G[^\x) + m^G^/\x) - g {[G[^\x)]'' + 3G[^\x)Gi^\x, x) + G^/\x, x, x)) = J{x). (2.15) 

We obtain the first Schwinger- Dyson equation for = 4 by setting J = 0: 

m^Gi -g[Gl + 3G'iG2(0) + ^3(0, 0)] = 0. (2.16) 

Using tlie same procedure as for = 3 we now take a functional derivative of Eq. ( p.l5| ) 
with respect to J{y) to obtain 

-d'G^'\x,y) + m'G^'\x,y) - g {3[G['\x)M'\x,y) 

+3G[^\x)Gi''\x,x,y) + 3Gi^\x,x)Gi^\x,y) +Gi^\x,x,x,y)) =6{x-y). (2.17) 

Now, setting J = gives the second of the Schwinger-Dyson equations for = 4: 

-d^G2{x -y) + m'G2{x - y) - g[3GlG2{x - y) 

+3G2(0)G2(a; - y) + 3GiG^{^, x - y) + 6*4(0, 0, x - y)] = 5{x - y). (2.18) 

Repeating this process once more by functionally differentiating Eq. (|2.17|) with respect 
to J{z) and setting J = gives 

—d'^Gsi^x — y,x — z) + rn^Gsi^x — y,x — z) 

—g[6GiG2{x — y)G2{x — z) + 3GlG3{x — y,x — z) 

+3G2{x - z)G3(0, x-y)+ 3G2{x - y)G3{0, x - z) + 3G2(0)G3(x -y,x-z) 
+3G'iG4(0, X - y,x - z) + G5{0,0,x - y,x - z)] = 0, (2.19) 

the third of the Schwinger-Dyson equations for N = 4. The fourth Schwinger-Dyson equa- 
tion is 

—d'^Gi{x — y, X — z, X — w) + rn^Gi^x — y, x — z, x — w) 

—g[6G2{x — y)G2{x — z)G2{x — w) + 6GiG2{x — y)G2,{x — z,x — w) 

+QGiG2{x — z)G3{x — y,x — w) + 6GiG2{x — w)G3{x — y,x — z) 

+3G'3(0, X — y)G3{x — z,x — w) + 3^3(0, x — z)G3{x — y,x — w) 

+3G3(0, X — w)Gs{x — y,x — z) + 3G^G'4(x — y,x — z,x — w) 

+36*2(3; — y)G4{0, X — z,x — w) + 3G2{x — z)G4{0, x — y,x — w) 

+36*2(2; — w)G4^{0, X — y,x — z) + 3G2(0)G4(x — y,x — z,x — w) 

+3G'iG'5(0, X — y,x — z,x — w) + Gq{0, 0,x — y,x — z,x — w)] = 0. (2.20) 

Again, we observe that the set of equations is incomplete; for this case the number of 
unknown Green's functions is two more than the number of equations (instead of one more 
as in the case A^ = 3). That is, the first equation for A^ = 4 contains the three unknowns 
Gi, G2, and G3, the second contains Gi, G2, G3, and G4, and so on. To solve this system 
of equations we truncate after the nth equation, but now to close the system of equations 
we must set Gn+i = and Gn+2 = in the last and next to last equations. 

How do we generalize this derivation to arbitrary A^? Clearly, when A^ is an integer and 
A^ > 2 we can use the two cases A^ = 3 and A^ = 4 discussed above as paradigms. All that 
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is needed is to functionally differentiate Eq. ( |2.7| ) N — 2 times. This gives an expression for 
(O|0^~"'^(a;)|O) in Eq. ( p.3|) . For example, when = 5 we have 



Z[J] {\gY\x)Y + Q[G^{'\x)fGi^\x,x)+^[G^^\x,x)f 

+AG[^\x)G[^\x,x,x) +G^^\x,x,x,x)) = (O|0^(x)|O), (2.21) 
when = 6, we have 

Z[J] {[G[''\x)]^ + 10[G[''\x)fGi''\x,x) + 15G[^\x)[Gi^\x,x)]'' 
+10 [G^-^^ (x)]^^^ ^ (x, X, x) + lOGi^^^ {x, x)Gi^^ {x, x, x) 

+5G[^\x)Gi^\x,x,x,x) +Gi^\x,x,x,x,x)) = (0|(/)^(a;)|0), (2.22) 

and so on. Once this calculation is completed, repeated functional differentiation with 
respect to J followed by setting J = gives the complete set of coupled Green's function 
equations. 

Note that while these equations are rather complicated, they are easily expressible in 
terms of multinomial coefficients. Following the notation of Abramowitz and Stegun [0, 
the multinomial coefficients are defined as follows: For each integer n, there is a set of 
multinomial coefficients; each coefficient expresses the number of possible ways to par- 
tition n = ai + 2a2 + . . . + na„ different objects into ak subsets containing k objects 
{k = 1, 2, . . . , n). For the first few integers n the sets of multinomial coefficients (called 
Ms in Ref. 0) are {1} for n = 1, {1,1} for n = 2, {1,3,1} for n = 3, {1,4,3,6,1} for 
n = 4, and {1, 5, 10, 10, 15, 10, 1} for n = 5. Observe that these are precisely the coefficients 
that appear in (^ for (O|0(x)|O), (p]) for (O|02(a;)|O), (CT ) for (O|03(a;)|O), for 



{0\(f)'^{x)\0), and (|2.22|) for {0\(f)^{x)\0) . Also, in these equations the powers of the Green's 



functions are the numbers ak- 

To be precise, for each n we must take all possible combinations of the numbers 
satisfying n = ai + 2a2 + . . . + nan- In the notation of Ref. vr = l°i, 2"^, . . . , n"". This 
allows us to read off the subscripts of the Green's functions and the powers to which they are 
raised; the numbers are the powers and the numbers they exponentiate are the subscripts 
of the Green's functions in that term. From this we calculate the quantity 

n' 

Mg = (n;ai,a2,...,a„) = — — " —- -, (2.23) 

(l!)''iai!(2!)"2a2! . . . (n!)''"an! 

which is the coefficient for that particular combination of Green's functions. 

For example, given (O|0^(x)|O) there are three possible combinations of the numbers a/, 
that satisfy n = ai + 2a2 + . . . + na„: ai = 3; ai = 1 and a2 = 1; and a^ = 1. This 
gives TT = 1^; TT = 1, 2; and vr = 3. Thus, there are terms of the form Gf, G1G2, and G3. 
Calculating M3 we find that the coefficients of these terms are 1, 3, and 1. Hence, we are 
able to reconstruct Eq. (|2.14|) . 



When N is noninteger, the situation is much more difficult. To derive an expression for 
(O|0^~^(x)|O) in this case, we assume first that is integer and use the general formula 
for (O|0^~^(x)|O) in terms of multinomial coefficients. Then, in principle, we can continue 
analytically off the integers using analytical expressions in terms of Gamma functions for 
these multinomial coefficients. 
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Of course, this procedure is complicated, but we illustrate it below by considering the 
simplest truncation possible in which we keep only the first two Schwinger- Dyson equations 
and set G„ = for all n > 2. In this case, for integer values of N we have the beautiful 
result that the only multinomial coefficients that appear are precisely the coefficients of the 
Hermite polynomials. Thus, for noninteger values of we have parabolic cylinder functions 
with the exponential divided off. It is in this fashion that we are able to continue off the 
integers n. 

To be explicit, we use the standard notation Dy{x) to represent the parabolic cylinder 
function Q and define the function Wy{x) by 

D^{x) = e-'''/^x''w^{x). (2.24) 

Then, for integer we factor [G^'^'']^ out of the equation for (O|0^(x)|O) and introduce the 
variables 'j^'^^ and 7 by 

^(j) ^ 1 ^ ^ and 7 = , ^ . (2.25) 



'Gi'\x,x) ' V^2(0) 

We then obtain for arbitrary noninteger 



00 



(This series terminates if is integer but is an infinite series if A^ is noninteger.) As we will 
see, when J = 0, the constant Gi is a negative imaginary number and the constant ^2(0) is 
real and positive. Thus, the argument 7 of is real and positive. 

Using Eq. ( p^.26D we write the field equation ( [2. 3] ) explicitly in terms of the Green's 
functions of the theory: 

- d^G^{'\x) + m^G[^\x) - gi''[G^{'\x)f~^WN-i{l^^^) = J{x). (2.27) 

(Remember that in the derivation of this equation we have discarded all Green's functions 
higher than ^2-) As before, we obtain the first of the Green's functions for arbitrary A^ by 
setting the source J = 0: 

m^Gi - gi^G^'^WN^i{l) = 0. (2.28) 

We obtain the second of the Schwinger-Dyson equations for arbitrary A^ by taking the 
functional derivative of Eq. ( |2.27| ) with respect to J{y): 

-d'Gi'\x,y)+m'Gi'\x,y) 

-^z^(^(Ar - l)[GS-^)(a:)]^-^G(-^)(x,y)i^;v-i(7^'^) 

+^[Gi^)(x)]^-^.^_,(7(^))^ffiML) = 5(x - y), (2.29) 

\/Gi'\x,x)J 
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where we have set ^3"^'' = 0. Now, setting J = in the above equation yields the second of 
the Schwinger-Dyson equations: 

—d^G2ix — y) + ni^G2{x — y) 



N 



G2 {x - y) 



5{x~y). (2.30) 



This equation can be simphfied by using the recurrence relations for parabolic cylinder 
functions [Bl, 



(2.31) 



to obtain 



- d'G2{x -y)+ m'G2{x - y) + (iV - \)g{iG^Y''w^.2{l)\G2{x - y) = 5{x - y). (2.32) 

Thus, Eq. (|2.28|) and Eq. (|2.32| ) are a closed system of two equations and two unknowns. 
The solution of this system for various choices of space-time dimension will be discussed in 
the following sections. 



III. ZERO-DIMENSIONAL THEORIES 



In zero-dimensional space-time the integral representation for the vacuum-vacuum func- 
tional Z[J] corresponding to the Lagrangian in Eq. ( |1.4| ) is an ordinary integral 



Z[J] 



exp 



2 ^ 



N 



+ J<p 



(3.1) 



To demonstrate the numerical accuracy of our truncation method we first study massless 
theories of the form in Eq. ( ^.1| ) where all quantities are exactly calculable. We find that for 
arbitrary the Green's functions can be expressed exactly in terms of Gamma functions. 
Then, we examine massive theories and show that weak-coupling diagrammatic methods are 
inadequate for the analysis of "PT-symmetric theories. 



A. Massless Theories in Zero Dimensions 



Using Eqs. ( |3.1|) , ( p.5|) , and (|2.6| ) the expectation value of the field (the one-point con- 
nected Green's function) is 



4ivy/"r(£+|)cos(f) 

9 J 



Gi = -i 

and the expectation value of the field squared is 
(O|02|O) fN 



(3.2) 



(0|0) 



9 J 



Hi) 




- 


-3cos2(^) 









(3.3) 
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Using Eq. ( p.8|) with J = we express the two-point connected Green's function as G2 = 
(O|0^(x)|O)/ (0|0) — G\, which can be calculated exactly using the two equations given above. 

Observe that for any > 2, Gi is pure imaginary and negative and G2 is real and 
positive, as claimed earlier. This is in fact true in any dimension. The reality of G2 and 
the imaginarity of Gi follows from the VT symmetry of these quantities. Applying V to Gi 
changes the sign, while applying V to G2 leaves the sign intact. Hence, under T, which acts 
as complex conjugation, Gi changes sign again and thus is pure imaginary, while G2 does 
not change sign and thus is real. 

The first truncation approximations to Gi and G2 obtained by keeping just the first two 
Schwinger-Dyson equations are found by solving the zero-dimensional massless versions of 
Eqs. ([2:2^ ) and {V^ . Observe that with m = 0, Eq. { ^J2E(\ demands that 



WN-i{l) = Q. (3.4) 

In zero-dimensional space-time there are no derivatives with respect to x and the delta 
function is unity. Thus, all Schwinger-Dyson equations are algebraic. The second Schwinger- 
Dyson equation is 

{N -l)g{iGr)''-^WM-2{lo)G2 = l, 

where 70 is the solution of Eq. (|3.4| ). (Note that there is some ambiguity with regard to 
which zero to choose. Our numerical studies suggest that the most positive zero is always 
the correct choice, but we do not have a proof. We find that it is this zero that gives the 
most accurate numerical results.) Recalling the definition of 7, we insert G2 = —G\l^1 to 
obtain an expression for G\. 



Gi 



{N - l)gwN-2{lo) 
Again using the definition of 7, we express G2 as 



l/N 

(3.5) 



G2 



(A^ - l)gwN-2{'yo) 



2/N 

(3.6) 



Table || compares the exact results with the corresponding first approximations for A^ = 3, 
A^ = 4, and A^ = 5. 

Deriving an expression for the second approximation in terms of A^ is difficult. However, 
for a specific A^, solving larger and larger systems of Schwinger-Dyson equations in which 
more and more Green's functions are included can be done symbolically on a computer. 
The case A^ = 3 is comparatively simple. If we keep n equations, then for k < n, the kth 
equation is linear in Gk+i- Thus, we can systematically solve a sequence of linear equations 
for the Green's functions, substitute the results into the nth equation, and solve the resulting 
polynomial equation for Gi. 

For A^ = 3, the first six approximations to iG^g^^^ are 0.79370, 0.69336, 0.74690, 0.71256, 
0.73987, and 0.71237. The exact answer, taken from Eq. (§3) with A^ = 3, is 0.72901 (see 
Table |). Note that the approximations oscillate around the correct answer. For the case 
of a Hermitian Hamiltonian we would expect that the approach to the correct answer is 
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monotone. This is because the calculation that we are performing is variational in nature. 
Keeping more and more Green's functions corresponds to enlarging the parameter space. 
However, in this case the Hamiltonian is not Hermitian, so this oscillation is not surprising. 
Indeed, it is consistent with our variational studies of PT-symmetric quantum mechanical 
systems in Ref. 0. 

Apart from the oscillations, we would hope that successive approximations would come 
closer to the correct value. Thus, the fact that the sixth approximation is worse than the 
fifth and that eventually most of the approximations become complex is a startling result. 
We believe that this is also due to the non-Hermiticity of the Hamiltonian. Indeed, while 
these approximations are numerically extremely accurate, we find that the convergence is 
rather slow and it may even be that this truncation method does not converge. To examine 
the convergence we have numerically solved very large systems of coupled Schwinger-Dyson 
equations for = 3; we have plotted in Fig. |l]all solutions for the one-point Green's function 
in the complex plane for the system of 150 coupled Schwinger-Dyson equations that lie near 
the exact answer. In this figure we see that the algebraic solutions for the one-point Green's 
function form a very small loop around the exact value. The largest dimension of the loop is 
along the real axis and is approximately 0.04. The furthest distance of a point on the loop 
from the correct answer is approximately 0.034. A weighted average of the points on the 
loop is extremely close to the correct answer 0.72901. 

The case A^ = 4 is more complicated than that of A^ = 3 because now the last two 
Schwinger-Dyson equations, rather than just the last equation, are nonlinear. As a result 
the largest number of equations we can solve and hence, the largest number of Green's func- 
tions we can include, is ten. The first five approximations to iGig^^^ are 1.10668, 0.90560, 
1.02988, 0.96159, and 1.02868. The exact value is 0.97774 (see Table |). Just as for A^ = 3, 
the approximations oscillate around the exact answer and successive approximations are nu- 
merically accurate but appear to converge very slowly. In this case, the fifth approximation 
is worse than the fourth. We have plotted in Fig. ^ all solutions of the first ten coupled 
Schwinger-Dyson equations for the one-point Green's function that lie near the exact an- 
swer. As in the case of Fig. Q these solutions form a small loop in the complex plane around 
the exact answer. 

The numerical work that we have done on zero- dimensional massless theories suggests 
that solving systems of truncated Schwinger-Dyson equations gives extremely accurate nu- 
merical approximations to the Green's functions. The convergence of the method is still not 
understood and warrants further study. 

B. Massive Theories in Zero Dimensions 

In zero-dimensional space-time, massless theories are so simple that it is possible to solve 
the Schwinger-Dyson equations as algebraic systems. However, in higher dimensions, the 
coupled Schwinger-Dyson equations are quite complicated. Thus, we will be interested in 
obtaining perturbative solutions for the case of small coupling constant g. 

To gain a rudimentary understanding of the perturbative nature of the theories in 
Eq. ( |1.4| ), we study in this subsection the path-integral expressions for the n-point dis- 
connected Green's functions F„ of a zero-dimensional —gcf)'^ field theory. These disconnected 
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Green's functions are expressed as ratios of integrals: 

_ JZo d<l> r exp (-|m202 + £04) 



dcj) exp |m202 + 



(3.7) 



Although these integrals cannot be done exactly, we can use the method of steepest 
descents to obtain an asymptotic result for small g with m fixed. We will see that, depending 
on whether the sign of m? is positive or negative, we obtain drastically different results. 

We begin by replacing by the dimensional quantity x according to = ^Ij g\m\x. 
This gives 

^ / 2H^\ r^da:a:"exp[-A(±x^-a;^)] 

\ g ) ro.rfa;exp[-A(±x2-x4)] ' ^"'•^^ 

where A = vtv" j g is a large positive parameter as ^ and we have distinguished between 
the two cases > and < 0. 

To calculate the integrals in Eq. (|3.7|) asymptotically, we find the stationary points of 
the function p{x) in the exponent, where 

p{x) = ±x^ - x^. (3.9) 

Setting the derivative of this function equal to zero gives the solutions a; = 0, ± 1/v^ for 
> and x = 0, ±i/^/2 for < 0. These stationary points are shown on Figs. |^ and ^. 
Next, we find the paths of steepest ascent and descent that pass through these stationary 

points. Substituting x = u + iv into Eq. ( ^.9[ ) and separating the real and imaginary parts 

yields 

p{u + iv) = ±{u^ - v^) - - t;^ + Gu'^v'^ + i[±2uv - ^uv[v? - v^)\. (3.10) 
Thus, the paths along which Imp = are given by -u = 0, f 7^ 0; f = 0, m 7^ 0; and 



u 



2 



±1/2, where the ± depends on the sign of rn? . 



To determine whether each of these paths is a steepest ascent or descent path, we perform 
a local analysis of the paths at each of the saddle points. The results of this analysis are 
displayed in Figs. |^ and ^. The stationary phase contour is chosen such that the boundary 



conditions are obeyed. Using Eq. (|1.3| ), we choose the end points that lie below the real axis 
at the angles — 7r/6 and — 57r/6. For > the contour follows the path u = — Jv^ + 1/2 



from the lower left quadrant of the complex-x plane up to the stationary point at — 1/a/2, 
then goes along the real axis from — 1/a/2 through the stationary points at the origin and 

at l/"\/2, and finally leaves the point 1/V2 and follows the path u = ^Jv'^ + 1/2 down to 
the lower right quadrant of the complex-x plane. For < the contour follows the path 
V = —\Ju'^ + 1/2, which passes through the stationary point at —i/\/2. 

Now, we determine the disconnected Green's functions F„ for both the m? > and 
< cases by evaluating the integrals in Eq. ( p.8| ) along the appropriate stationary phase 
contours. For > the Green's functions with odd subscript vanish (by oddness) along 
the part of the contour that lies on the real axis. Evaluating the integrals asymptotically 
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along the remaining parts of the contour gives an exponentially small result; apart from a 
multiplicative constant 

F,^^,^g—'l\-'^/^ (g^O^). (3.11) 

The dominant contribution to Green's functions of even subscript comes from the portion of 
the contour along the real axis. Evaluating the integrals near the saddle point at the origin, 
we find that the small-^f behavior of these Green's functions is given by 

F2n~m-2" (^?^0+), (3.12) 

apart from a multiplicative constant. Note that the Green's functions with even subscript 
behave very differently than the Green's functions with odd subscript. Further, notice that 
the behavior of the Green's functions with odd subscript is inherently non-perturbative 
and, as a result, Feynman perturbative methods cannot be used to calculate these Green's 
functions. 

For < the dominant contribution to all of the Green's functions comes from the 
saddle point at — i/y^. As a result, all of the Green's functions exhibit similar behavior. 
To be specific, we have 

Fr,r^{~ir\mrg-^/' ig^O+). (3.13) 

Thus, once again, it is clear that these results cannot be obtained using Feynman diagrams. 
As we will see in Sec. |V| this theory is also the one that is asymptotically free in four 
dimensions and the more interesting of the two cases. 



IV. NUMERICAL STUDY OF ONE-DIMENSIONAL THEORIES 

Our success with zero-dimensional massless theories prompts us to study one-dimensional 
(quantum mechanical) massless theories of the form in Eq. ( |1.4|) . Numerical computations 
have been performed in one-dimension which allow us to compare with exact numbers |]l|. 

In analogy with the last section, we find the first approximation to Gi and G2 by solving 
the one-dimensional massless versions of Eqs. (^^281) and (lO^) . Observe that with m = 



we obtain Eq. ( p.4| ) once again. In fact, Eq. ( |3.4| ) holds independent of the dimension. The 
second Schwinger-Dyson equation is given by 

- d'G2{x -y) + {N- l)g{iGif-'wN-2{lo)G2{x - y) = 5{x - y). (4.1) 

This equation depends on the dimension D through the partial derivative. If we introduce 
the variable M defined by = {N -l)g{iGi)^-'^WN-2{lo), it is clear that Eq. (|1]) is just 
the equation for the Feynman propagator, whose solution in one-dimension can be written 

Consequently, 6*2(0) = 1/(2M). Recalling the definition of 7 we use 6*2(0) = —G\/'^'^ as we 
did in the previous section to obtain an expression for Gy. 
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_A{N - l)gwN-2ilo) 
Substituting Gi into our expression for M yields 



l/iN+2) 



(4.3) 



M= {N -l)gwN-2{lo) 
which further allows us to write 6*2(0) as 



Af-2 



2/iN+2) 



(4.4) 



^2(0) = - 



I / 70 



2-N 



2/iN+2) 



(4.5) 



For field theories, M represents the renormalized mass, which is nothing more than the 
difference in energy between the first excited state and the ground state Ei — Eq. Table |I| 
compares the exact values of M and Gi with the corresponding first approximations for the 
cases = 3, = 4, and = 5. (Here, we must set g = N/2 and multiply M by 2 to 
match the Lagrangians studied in Ref. IQ.) 

In one dimension it is difficult to obtain the second approximation, even for a specific 
N, because it requires solving coupled systems of nonlinear differential equations. As in 
the previous section, the easiest case to study is = 3. The first three Schwinger- Dyson 
equations are given in Sec. | as Eqs. (|2.10|) , (p.l2|) , and (|2.13| ) with m = 0. To close this 
system of equations we set G^ = and z = y. Then, the third Schwinger- Dyson equation 
becomes 

- d^Gsix -y) + gt[2GiGs{x - y) + 2Gl{x - y)] = 0. (4.6) 
Next, we Fourier transform the second and third Schwinger- Dyson equations to obtain 



ip' + M')G2ip) + giG.ip) 



and 



{p^ + M')Gs{p) + 2gz 



^ fin ~ 

7rG2{q)G2{p - q) = 0, 
00 Zn 



(4.7) 



(4.8) 



where we have used the same definition of M as above, the convolution property of Fourier 
transforms, and the translation invariance of Green's functions (^3(0, x — y) = G^{x — y^x — 
y) =G^{x-y). 

We now solve for G2{p) and obtain 



G2{P) 



1 



2g' 



p2 + M2 (p2 + M2)2 



30 rjn _ 

^G2{q)G,{p-q). 

00 ZTT 



(4.9) 



The simplest approach to solving this nonstandard integral equation is to iterate it to 
high order in the small parameter g. This iterative procedure can be represented in terms 
of diagrams. These diagrams all have a similar structure: They begin with one line that 
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branches into two lines. This branching process continues until the maximum number of 
lines is attained. Then the process is reversed, with lines combining in pairs until only one 
line remains. We were able to perform the calculations symbolically on a computer. We 
calculated the propagator to order 0(5'^°); the first three terms in the expansion are 

G2(P) 



(456M^ + 70M2p2 + V)^4 

9M6(p2 + M2)2(p2 + 4M2)2(p2 + 9]Vf2) ^ ^ - ) 



+ 



Now that we have 6*2 (p) to high order in we make the ansate that it can be expressed 
in the form 

^^(^) = ^^ + ^^ + ^^ + ---' ^'-''^ 

where M„ = nM + hi^ng"^ + h2,ng^ + &3,nfi'^ + • • ., Zi = 1 + ai^ig'^ + a2,ig^ + 03,15'^ + • • •, 
Z2 = ai,25'^+o.2,2fi'^+fl3,25'^ + - • •, ^3 = ^2,35'^+ 03,35'^ + • • •, and so on. By matching this ansatz 
to our calculation, we determine the coefficients ak^n and bk^n- [The expansion of the ansatz 
does not exactly match the expansion of our calculation. This is easily understood because 
our calculation for G2{p) only involves special diagrams described above, while the ansatz 
involves all types of diagrams. However, the system of equations is neither overdetermined 
nor underdetermined, so all coefficients may still be calculated.] 
The series for the M„ are 

_ g^ _ 31^^ 1279/ 98287/ 9641179/° 

^~ ^ 3M^ ~ 72M9 ^ 1944M14 ~ 93312M19 ^ 5598720M24' 

^ 3M4 108M9 1944M14 279936M19 ' 

216M9 486M14' 

4/ 109/ 

M4 = 4M + + 

3M4 270M9 ' 

5o2 

M^ = 5M + -^. (4.12) 

^ 3M4 ^ ' 

Since we now have expressions for M„ in terms of M, we need only determine M ac- 
curately to finish the calculation. This is done using the first Schwinger-Dyson equation 
( p.lOp , which implies that 

= 4/^2(0). (4.13) 

Observe that 

^2(0) = J^^^G2{p). (4.14) 
So, based on our calculation of the first three terms above, the first few terms are 
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1 q2 74 

^2(0) = — ^ + (4.15) 



This allows us to express M and (7 as 




' +wA7TT + 0{9'). (4.16) 



9M6 24M11 

Keeping terms to order 0{g^^) we obtain M = 1.126151(7^/^. 

With this value of M, the M„ become numerical series multiplied by an overall factor 
of g'^^^. Successive terms in this numerical series decrease in magnitude and thus the series 
appears to be convergent. To compare with the numerical results and match the Lagrangian 
in Ref. we set (7 = 3/2 and multiply the Mn by 2; results are given in Table |ITJ 



V. SCHWINGER-DYSON EQUATIONS IN D DIMENSIONS 

In this section we show how to solve truncated systems of Schwinger-Dyson equations in 
D dimensions. First, we consider the case D < 2, in which it is not necessary to perform 
any renormalization. Then we consider the case of arbitrary D, in which it is necessary to 
discuss renormalization. 



A. Schwinger-Dyson Equations for D < 2 

In this subsection we solve Eqs. (|2.28|) and ( |2.32| ) in arbitrary dimension D with m = 0. 
This calculation is a straightforward generalization of the one for D = 1. 

As previously stated, when m = 0, Eq. ( |2.28| ) implies Eq. (|3.4|). In addition Eq. ( p.32| ) 
continues to imply Eq. 



- d^G2{x -y)+ M^G^ix -y)= 6{x - y), (5.1) 

in which we have defined the renormalized mass by 

M' = -gi^'iN - l)G^-'wM-2{lo). (5.2) 

Equation (|5.1|) is just the differential equation satisfied by the Feynman propagator. 

We solve these two equations in arbitrary dimension D by taking the Fourier transform 
to obtain the propagator: 

G2(p) = l/(p2 + M2). (5.3) 
Fourier transforming this propagator back to position space and then setting x = y gives 

^2(0) = M^-2r(l - D/2)7r-^/22-D_ ^5 4) 

Using Eqs. ( |2.25| ) and ( |5.2D we solve for Gi. 
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= -^{[iN - l)gw^-.{io)r''y'ril - Z^/2)(4.)-^/^o^}^/^-™™^ . (5.5) 
Substituting this result into Eq. ( |5l^ ) yields 



and substituting the last two results into Eq. ( ^.41 ) gives an expression for G2- These expres- 
sions reduce to the zero-dimensional and one-dimensional solutions given in the previous 
two sections. Notice that each of these expressions becomes singular at D = 2. 

Using these approximate solutions we can determine the large- behavior of the Greens 
functions as ^ oo. To do so, we need to determine the asymptotic behavior of the 
zeros of the parabolic cylinder function. According to Ref. 0, the largest zero of Wn-iij) 
for large n is given by 70 ~ '2^/n. Substituting this into the integral representation for 
u!n-2{'y) and using Eq. ( |2.31| ) we perform a steepest-descent calculation to obtain Wn-2 ~ 
y^2/7r3~"^/^r(l/3)2~"'n"^''^e"/^. We have verified these results numerically to high accuracy. 
Taking A^ — > 00 in the expressions for Gi, G2, and M above and using the asymptotic results 
for the parabolic cylinder function, we obtain 



Gi 2i/Ve = -1.21306i, 

G2 ~ l/{Ne), 



er(l-D/2)(47r)-^/2^l'^^' (5.7) 



Observe that Gi is independent of both A^ and D to this order, G2 depends on A^ but 
not on D, and depends on both A^ and D. These results are valid for large A^ for all 
D < 2. These properties are evident in zero dimensions from Eqs. (|3.2| ) and ( p.3|) . While the 
behavior of Gi and (O|0^(O)|O) is correctly predicted, the behavior of G2 is not because the 
first-order behaviors of Gi and (O|0^(O)|O) cancel, leaving a second-order term to describe 
6*2 • Also, this calculation predicts that M increases like Ar2/(2-D)^ r^Yiat is, for large A^ the 
separation between the energy levels diverges and hence, the energy levels must diverge. In 
D = 1 our approximation suggests that M grows like A^. In fact, M grows like A^^ in D = 1 
as discussed in Ref. [0. This is the simplest possible truncation, but it suggests the correct 
behavior. 



B. Perturbative Renormalization of the Schwinger-Dyson Equations in 
/^-Dimensions and Leading-Order Calculation of the Beta Function 

Let us first consider the "PT-symmetric ig(j)^ theory in six space-time dimensions. The 
Green's functions for that theory are governed by the system of equations beginning with 
Eq. (|2.10|) , ( p.l2| ), and (|2.13|) . Let us seek a perturbative solution to this system of equations 
in which Gn ~ g"'~'^- In leading order we have first, from Eq. ( p.lOj ), 

2 

Gi = t—. (5.8) 

g 

Then, from Eq. ( |2.12| ), the two-point function in momentum-space is 
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where 



= m^ + 2igGi = -m^ 



1 



p2 + M2 ' 



(5.9) 



(5.10) 



Thus, in order to avoid unphysical singularities, we must have for this type of solution 



< 0. The leading solution to Eq. ( |2.13D is 



2ig 



(p2 + M2)(g2 + M2)[(p + g)2 + M^] ' 



(5.11) 



which has an obvious interpretation as a vertex with three external lines. 
More generally, the solution to Eq. (|2.1(]| ) is 



Gi 



- yim'/g ± i^{m^/gf + AG2{Si)\ , 
which corresponds in the perturbative case, where 1(72(0)1 <^ [rn^/gY, to 

Gl = { 3 ' 



(5.12) 



(5.13) 



The second solution given in Eq. ( p.l3| ) corresponds to the usual perturbative tadpole con- 
tribution to the vacuum expectation value of the field, while the first is the new, nontrivial 
solution given in Eq. (|5.8|). 

Perturbatively solving the next in the sequence of Schwinger-Dyson equations, we obtain 
for the four-point function in leading order 



Gi{p,q,r) 



-2ig 



{{p + q + r)2 + M2)(p2 + M2)(g2 + M2)(r2 + M2) 
1 1 1 



X 



+ 



+ 



(p + g)2 + M2 (p + r)2 + M2 (g + r)2 + M2 



(5.14) 



Inserting this back into Eq. ( |2.13|) , we obtain the one- loop correction to the three-point 
function. Apart from a tadpole term, this is just the same as that found in the conventional 
0^ theory; correspondingly, the beta function is obtained from that in the conventional 
theory by the replacement g ^ ig (see, for example, Ref. WW l 



(5.15) 



where g{n) is the running coupling at scale /i. Unlike the usual 0g theory, the VT symmetric 
theory given here is not asymptotically free. 

Of course, the (pf theory is of far greater interest. In particular, it plays a crucial role in 
the standard model as the origin of particle masses through the Higgs mechanism. Yet the 
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triviality of that theory is a source of difficulty. What happens here, when we set = 4 in 
Eq. 



The first few Schwinger-Dyson equations are given in Eqs. ( |2.16| ), ( p.l8| ), ( |2.19| ), and 
( p.20| ). Note that the last three equations can be simplified through the introduction of the 
renormalized mass 

= - 3gG2{0) - 3gGl (5.16) 

Thus we obtain the following equations for the two-point function, 

i-d' + M^)G2{x -y)- gGiG,{Q, x-y)- (7^4(0, 0, a; - y) = 5{x - y), (5.17) 

the three-point function 

(-9^ + M^)G3{x -y,x-z)- QgGiG2{x - y)G2{x - z) - 'ig{G2{x - i/)G3(0, x - z) 
^G2{x - z)G'3(0, x-y)\- 3(7G'iG4(0, x-y,x-z)- gG^{^, 0, x - y, a; - z) = 0, (5.18) 

and the four-point function 

(—(9^ + M'^^G/^{x — y,x — z,x — vS) — 6gG2{x — y)G2{x — z)G2{x — w) 

— 3g[G2{x — y)G4{0, x — z,x — w) + G2{x — z)G4^{0, x — y,x — w) 
+ G2{x — w)G4{0, X — y,x — z)] — 6gGi[G2{x — y)G^{x — z,x — w) 
+ G2{x — z)Gs{x — y,x — w) + G2{x — w)G3{x — y,x — z)] 

— 3g[G3{0, X — y)G^{x — z,x — w) + 6*3(0, x — z)G^{x — y,x — w) 
+ 6*3(0, X — w)Gs{x — y,x — z)] — 3gGiG5{0, x — y,x — z,x — w) 

— SgGei^O, 0,x — y,x — z,x — w) = 0. (5.19) 

Now there are two regimes. If > the only consistent perturbative solution to the 
above system of equations is one in which the odd Green's functions are exponentially small, 
6*271+1 ~ e^^/^, and the even and odd Green's functions decouple. [This result is analogous 
to that in Eq. ( p.ll|) .] The even Green's functions possess the same perturbative expansion 
as in the usual 0^ theory except for a change of sign of the coupling constant, so again the 
sign of the beta function reverses: 

g^2 



(^(9) = -27[^^) . (5.20) 

This theory is asymptotically free. Furthermore, in the nonperturbative regime it exhibits 
parity symmetry breaking, but possesses VT symmetry because Gi is imaginary. 

The other regime is even more interesting. If < it is consistent to proceed in analogy 
with our treatment of the (j)^ theory above. We may assume a perturbative solution of the 
form Gn ~ (7"/^^^, as we have already seen in Eq. (p.l3|) . Then we have a purely imaginary 
vacuum expectation value of the field, from Eq. ( |2.16| ): 

Gi = (5.21) 
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while the leading two-point function has the usual form of a propagator: 



= (5.22) 

Here the renormalized mass in leading order is positive: 

= m^ - 3gGl = -2ml (5.23) 
The leading three-point function has an evident diagrammatic interpretation: 

= '^ %^ + M^g^ + MMp + g)'+M^ - ^'-''^ 
The tree-level four-point function is easily extracted from Eq. ( |5.19| ): 

6^ 



3m^ 3m^ 3m^ 



1+ . . 



(5.25) 



being composed of contibutions from primitive four-point and three-point vertices. 

Now we have perturbative parity symmetry breaking: the scalar field acquires a vacuum 
expectation value comparable to that of the gauge bosons in the standard model, y^Gi. 
Further, it appears likely that the theory is asymptotically free, because the sign of the 
four-point vertex is reversed. Indeed, apart from one-particle-reducible graphs, Eq. ( ^.25] ) 
gives just the usual primitive vertex in the high momentum limit, except for a change in 
sign. The theory is renormalizable because, apart from divergences associated with the 2-, 
3-, and 4-point functions, no additional divergences occur. This is due to the fact that, for 
example, the 5-point function has no primitive vertices, as can be easily seen from the next 
in the sequence of Schwinger-Dyson equations after Eq. ( [5.181 ). The lowest order diagrams 
contributing to G5 are as sketched in Fig. ||. 

We have shown that the signs of the beta functions for a conventional gcf^ theory and for a 
PT-symmetric igcl)^ theory in six space-time dimensions are reversed. Thus, while the former 
theory is asymptotically free, the latter is not. Similarly, the beta functions for a conventional 
gcj)'^ theory and for a PT-symmetric —gcf)^ theory in four space-time dimensions are reversed. 
Thus, while the former theory is not asymptotically free, the latter is. Similarly, as we have 
already argued in Ref. |T2|, we believe that while conventional quantum electrodynamics is 



not asymptotically free, PT-symmetric electrodynamics is asymptotically free and possesses 
a nontrivial fixed point. 

One of us, CMB, thanks Arthur Lue for many useful discussions at the Aspen Center 
for Physics. This work was supported in part by the U.S. Department of Energy. 
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FIG. 1. Solutions of the first 150 coupled Schwinger-Dyson equations for the dimensionless 
one-point Green's function iGig^^'^ for the case of a massless zcro-dimensional = 3 theory. Note 
that the solutions (indicated by dots) lie in a small portion of the complex plane very close to the 
exact answer 0.72901. 
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FIG. 2. Solutions of the first ten coupled Schwinger-Dyson equations for the dimensionless 
one-point Green's function iGig^^^ for the case of a massless zero-dimcnsional = 4 theory. The 
solutions (indicated by dots) lie in a small portion of the complex plane very close to the exact 
answer 0.97774. 
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FIG. 5. Lowest order graphs contributing to G5. 



24 



TABLES 



TABLE L Exact values iC^^^^g^l^ and Gf^^g^/^ [see Eqs. (|]|) and (|3j)] compared 
with the first approximations iOf^g^^^ and G^^g'^^^, which are obtained from the first two 
Schwinger-Dyson equations and given in Eqs. ( |3.5| ) and (|3.6D. This is done for the three cases 
A'^ = 3, = 4, and N = 5. Observe that the percent error increases with A^ for Gi and decreases 
with A^ for G2. These trends continue until Gi has a maximum error of 23.2% at A^ = 53, and 
until G2 has a minimum error of 0.34% at A^ = 8. For A^ > 53 the error for Gi decreases until 
it levels off at 21.3%. For A^ > 8 the error for G2 increases like A^. The large-A^ behavior of our 
approximations is more fully discussed later. 
N iGf^'g^l^ igSPgi/jv 

3 0.72901 0.79370 0.53146 0.27516 

4 0.97774 1.10668 0.28000 0.14907 

5 1.07865 1.24829 0.16433 0.10158 



TABLE IL Exact values of iGf^* and W'^'^^ = Ei-Eq (see Ref. |l|) compared with the first 
approximations iG^ and M^^, which are obtained from the first two Schwinger-Dyson equations 
(|4.3| ) and ( |4.4| ). This is done for the three cases A^ = 3, A^ = 4, and N = 5. Note that the error is 
worse than for the zero-dimensional theories, and that the error for M increases with A^. However, 
the error for iGi is smallest for A^ = 4. 

A^ iGf^^ iGf^ ^exact j^jSD 

3 0.59007 0.37011 2.95293 2.70192 

4 0.86686 0.82548 4.52620 3.63424 

5 1.01310 1.15416 6.70000 4.72160 



TABLE III. 

Schwinger-Dyson approximations for M„ compared with exact values of the energy differ- 
ence, En — Eo, calculated in Ref. |jl[. These approximations are based on the truncation of 
the first three Schwinger-Dyson equations for an ig(f)^ field theory of the type in Eq. (1.4). 



Notice that Mi is greater than the numerical result in this case, while in Table ||, Mi was 
less than the numerical answer. This suggests that in one dimension the oscillatory nature of 
successive approximations is present once again. Moreover, the percent error has decreased 
significantly; the error is 8.50% for Mi calculated using the first two Schwinger-Dyson equa- 
tions compared with 3.52% calculated using the first three Schwinger-Dyson equations. 



n 



En — El 







M, 



SD 



/o error 



1 
2 
3 
4 
5 



2.952962 
6.406007 
10.158155 
14.135286 
18.295263 



3.056763 
6.191828 
9.610135 
12.871434 
15.681836 



3.52 
3.34 
5.39 
8.94 
14.28 
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